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We have studied the thermal properties of atom chips consisting of high thermal conductivity Aluminum 
Nitride (AIN) substrates on which gold microwires are directly deposited. We have measured the heating of 
wires of several widths and with different thermal couplings to the copper mount holding the chip. The results 
are in good agreement with a theoretical model where the copper mount is treated as a heat sink and the thermal 
interface resistance between the wire and the substrate is vanishing. We give analytical formulas describing the 
different transient heating regimes and the steady state. We identify criteria to optimize the design of a chip as 
well as the maximal currents Ic that can be fed in the wires. For a 600 /xm thick-chip glued on a copper block 
with Epotek H77, we find /c = 16 A for a 3 /xm high, 200 /xm wide-wire. 

PACS numbers: 39.25.+k, 03.75.Be 



I. INTRODUCTION 

A few years after the first Bose Einstein condensates (BEC) 
in atomic vapours, the will to miniaturize the set-ups has led 
to the production of atom chips, which consist in microfabri- 
cated elements, most usually wires, that are used to trap and 
manipulate cold atomic clouds. The atom chip soon became 
increasingly popular as a compact, robust and versatile device 
suitable for the production of BEC 111 Q, [sl] and for studies 
on quantum matter, cold atom-based metrology, or quantum 
information |0] . 

The major interest of atom chips is to manipulate atoms 
in the close vicinity of the field sources, so that trapping po- 
tentials with strong spatial variations can be obtained at very 
low power consumption. Small structures and tightly con- 
fining traps allow one to perform efficient evaporative cool- 
ing at high collision rates, to squeeze cold clouds to very 
anisotropic geometries so that low dimensional regimes are 
reached d, [g, S S], and make it possible to realize a great 
diversity of trapping geometries 1 4] . 

Still, there are limits to the miniaturization of the structures 
and for several reasons, large currents can still be needed. 
First, as the loading stages require traps deep and wide enough 
to collect a high number of atoms 13], big structures running 
large currents are needed. These large loading structures can 
be placed below the chip ifToll , but other limitations still pre- 
vent the use of arbitrarily small structures. One problem is 
that when atoms are brought close to the source wire they 
become sensitive to the potential roughness created by the 
wire's imperfections (Til [12]. Second, it is often desirable 
to avoid the interactions between the atoms and the surface 
(Casimir-Polder force (illlH] or Johnson noise causing spin- 
flip losses [IB]), which involves keeping the atoms at least 
some microns away from it. For all those reasons, the maxi- 
mal current that can be carried by each wire is a crucial param- 
eter determining the possibilities of an atom chip. Unless one 
uses superconducting wires 1 17], in which case the maximal 
current is the critical current above which the metal becomes 
normal, the maximal current is determined by the dissipation 
of the heat generated in the wires. In all this paper, we only 
consider resistive wires, much simpler to achieve experimen- 
tally and most widely used. 



Up to now, atom chips have mainly been realized on Sil- 
icon (Si) wafers because Si is cheap, fabrication techniques 
are well developed, and it is a good thermal conductor. But, 
as Si is semiconductor, an electrically insulating layer, gen- 
erally Si02, needs to be placed between the wafer and the 
metallic wires. Unfortunatel y, S i02 is also a thermal insula- 
tor, and, in a previous study llisll . it was found that this layer 
is the main limitation to the removal of heat in Si-based atom 
chips. On the other hand, AIN is a substrate material that 
has been especially selected for being simultaneously a good 
electrical insulator and a good thermal conductor. It was first 
used for high-power microelectronics applications |[l9l] . More 
and more groups working on atom chips are now moving to 
AIN substrates H HO, [13 [H, because they allow for direct 
deposition of the wires on the substrate. Since no thermal con- 
tact resistance between the wire and the substrate is expected, 
much better heat dissipation is foreseen. 

We have fabricated such chips and measured their thermal 
behavior. We first present the model that we have developed to 
understand our experiments and identify the different heating 
regimes. We then show that this model reproduces very well 
the measured wire heating for different wire widths, different 
thermal couplings to the copper mount and within the differ- 
ent heating regimes. In particular, we show that the heating 
in AlN-based atom chips is only governed by heat diffusion 
in the substrate, unlike the Si-based atom chips. Therefore, 
the relevant phenomenon to consider for the current limita- 
tions is the long-time heating rather than a fast heating due to 
a thermal resistance between the wire and the substrate. The 
thermal coupling to the copper mount and the possibility to 
quickly reach a stationary regime thus acquire a crucial im- 
portance. Finally, we draw pratical conclusions on the optimal 
design and the current limitations of AlN-based atom chips. 



II. THEORETICAL MODEL 

At any time of its operation, the temperature of a current- 
carrying wire is the result of its resistive heating and the heat 
removal via different channels. Most of the heat is removed 
by conduction through the substrate on which the wire is de- 
posited, out to the chip mount, considered as a heat reservoir. 
When the chip is placed in a vacuum chamber, the only other 
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A 


AIN conductivity 


128 W/(Km) 


D 


A/c : AIN diffusion constant 


53 mm^/s 


a 


thermal contact resistance between 
substrate and heat sink 


1.1 X 10"^ m^KAV 


e 


substrate thickness 


600 fim 


Ic 


Act : contact length 


14 mm 


Is 


Veh : stationary length 


2.9 mm 




/ D — crec '. stationary time 


0.16 s 


te 


/ D : crossover time between 2D 
and IDx regimes 


6.7 ms 




resistive power generated by the wire 
per unit length 


920 Wm-^ 



TABLE I: Relevant parameters for the wire heating and their value in 
the case experimentally studied in this paper where the chip is glued 
onto the copper block. The value of $ is given for a 200 /xm wide, 
3 fim high gold wire running a current of 5 A. 



mechanism for heat dissipation is black-body radiation of the 
wire which, as we show later, has a negligible effect. In air, 
as for some of the experiments presented below, air convec- 
tion can also play a role, but it is expected to be negligible, in 
particular since the wires in the chip we use are covered by a 
6 /im thick-layer of resist of high thermal resistance. In the 
following, we thus assume that heat conduction inside the sub- 
strate is the only mechanism for heat removal and we compute 
the expected heating of the wire. 



Au wire 




FIG. 1 : Sketch of our model. The substrate is infinite in the x direc- 
tion and the Cu mount acts as a heat reservoir. The heat dissipated 
in the wire diffuses in the substrate in the IDz, 2D and IDx regimes 
successively. When the transverse spread is h, an equilibrium state 
is reached (see text), (color online) 



In our model, depicted in Fig. [TJ we assume that the wire, 
of width W along the x direction is infinite in the y direction. 
We also assume that the substrate of thickness e is infinite in 
the transverse direction x. The back surface of the substrate is 
supposed to be in contact with a heat reservoir at temperature 
To with a thermal contact resistance a (in m^KAV). Finally, 
we assume that the energy flux per unit area from the wire 
to the substrate is homogeneous over the wire width. Using 
these hypothesis, we compute the temperature inside the sub- 
strate and in particular just below the wire. We assume that 
the thermal contact resistance between the wire and the sub- 
strate is vanishing so that the substrate temperature at the wire 



position is equal to the wire temperature. 

The substrate thermal properties are described by its spe- 
cific heat per unit volume c and its thermal conductivity A. 
The resulting diffusion constant is D = A/c. With the thermal 
contact resistance between the substrate and the heat reser- 
voir, we can construct a length Ic = Xcr. However, a more 
relevant length scale is Ig = v^el^, denoted as the stationary 
length, as we explain in the following. Thus, three different 
length scales govern the heat diffusion process : W, e and Is. 
From these lengths we derive the corresponding time scales 
tw = W^/D, te = e^/D andt, = I'^JD. Table U gath- 
ers the most relevant parameters and gives their value for the 
experimental case studied in this paper, in the most favorable 
situation where the chip is glued onto the copper block. 

Let us first consider the step response of the system for a 
wire of constant resistivity : the wire and substrate are at tem- 
perature To when the current is suddenly turned on, producing 
a heat flow from the wire to the substrate per unit length and 
per unit time <l>. We compute the evolution of the wire temper- 
ature, noting by AT its deviation from To. In the equations, 
we express temperature in energy unit, setting ks = 1. To 
give some physical insight we consider the situation where 
the characteristic times satisfy tw ^ te <C t^, a situation 
usually fulfilled. For times t <^ tw, the energy has diffused 
inside the substrate over lengths much smaller than the wire 
width W. Then the diffusion inside the substrate is expected 
to be one-dimensional in the z direction (see Fig.[T]) and in this 
IDz regime the wire heating can be written as 
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Next, for times tw ^ t <^ t^ the energy has diffused into 
the substrate over lengths much larger than the wire width but 
much smaller than the substrate thickness e, so that the heat 
diffusion is well described by a two-dimensional model (see 
Fig.[T]) and we expect the temperature to increase as 
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where VFeff is an effective width on the order ofW. At times 
t ^ te, the energy has diffused over lengths much larger 
than the substrate thickness e, so one expects now a one- 
dimensional model in the x direction (see Fig. [T]) to be an 
accurate description of the wire heating, which takes the form 
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This formula is valid as long as the energy flux per unit length 
to the reservoir ^r is much smaller than The heat flux to 
the reservoir per unit area is jr = {T{x)—To)/a where T(x) 
is the substrate temperature at position x. Using T(0) — To = 
AT and since T{x) — Tq falls off to zero over a width of order 
L Voi, we find that = / dxjR AT^/Dt/a. Insert- 
ing in Eq.El we find ^r ^ t^D/{elc) = ^t/ts. Thus Eq.[3] 
is valid as long as t <C ts. For t ^ tg, ^r ^ : the heat 
flux to the reservoir compensates the input heat flux <l> and the 
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FIG. 2: Transient wire heating for = 3 x 10~'^/c and e = 3 x 
10~^/c (dashed hne). Temperature and time are scaled to $/A and 



tc 



r. I D respectively. The solid lines correspond to fits with the 
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functions (from earlier to later times) ai * AT ^, AT 
bsAT^^"", the fitted parameters being ai = 0.74, Vl4ff 
as 1.0 and 0.79. 



and as + 
: 0.17W, 



librium temperature is then obtained setting to in Eq. [S] 
In the case Ic ^ e W, the wavevectors much smaller than 
1/e give a lorentzian contribution, whose integral yields the 
prefactor a = 1 /2 for the first term in rhs of Eq. |4l 

Most importantly, we use Eq. \5\ to compute the transient 
heating of the wire. Fig. [21 shows the result of the calculation 
for an input energy flux ^ which is a step function in time, 
for well separated length scales. Fits in the three different 
temporal domains are in good agreement with the analytical 
laws derived above. We observe that the initial IDz behavior 
fails noticeably already for t ^ and that Wgff in Eq.[2l 

is noticeably smaller than W (Weff - 0.17(3)W^). On the 
other hand, the transition from the 2D to the IDx regime and 
the convergence to the equilibrium state occur as expected for 
t c^te and t c^ts respectively. 

The above calculations assume that the energy flow <l> is 
constant in time, however, as the wire heats up, its resistivity, 
and consequently <l> as well, also increase. To take this effect 
into account numerically, we compute the wire temperature as 



system reaches a steady state. This justifies the denotation of 
ts and Is as the stationary time and length respectively (see 
Fig.[T]). Interestingly, tg, which can be rewritten as tg = cea, 
does not depend on A. Note also that tg is equal to the re- 
laxation time of the substrate, defined as the time constant of 
the exponential decay of the substrate temperature after a ho- 
mogeneous heating. This can be understood by the intuitive 
picture that the heat spreads in the substrate as long as it has 
not been absorbed by the mount. 

When the stationary flow is finally established, the equilib- 
rium temperature at the wire can be estimated by summing up 
the contributions of the three successive regimes. The con- 
tribution of the late IDx regime is found replacing t = tg 

in Eq. [3l T^^^^ = Q^x^^, where a is a numerical prefac- 
tor. From Eq. [21 we find that the 2D regime contributes to 
an amount T^^ = where Wq is of the order of 

W. As for the initial IDz regime, its contribution scales as y 
so that it is possible to incorporate it in T^^ by renormalizing 
Wo . Finally we get the analytic expression for the equilibrium 
temperature 
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where we have replaced a by ^ a posteriori from the calcula- 
tion that we explain below. 

In order to get quantitative predictions and to investigate 
the crossovers between the different diffusion regimes, it is 
useful to perform numerical calculations. For this purpose we 
use standard Fourier analysis. For an input heat flow at the 
top of the substrate that varies as jfe,a;e*^^e*^^ we find the 
temperature at z = 



AT,. 



'■ Jk,c 



Sinh(Ke) /Ic + KCosh(Ke) 1 
' Cosh(Ke)/lc + KSinh(Ke) Zx' 



(5) 



where K = Vk'^ -\- icXuo. We make the assumption that the 
heat flow j(x) is homogeneous over the wire size. The equi- 



T(t) 



R(t-t')^(t')dt' 
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where the impulse response function R{t') is determined us- 
ing Eq.[5l To compute ^{t') at a given t^ we assume a homo- 
geneous wire temperature, so that ^{t') = P p{T{t')) /hW . 
This approximation is justified a posteriori since the computed 
wire temperature inhomogeneity never exceeds 10% for the 
parameters explored in this paper. The dependence of the gold 
resistivity with temperature is approximated by 



/) = po(l + c^AT) 



(7) 



where a = 0.0038 is obtained from a linear fit of the 
reported values of p between 200 K and 500 K ll23ll . 

Our last remark about the model concerns the assumption 
that the input energy flux j{x) is homogeneous over the wire 
width. This is expected to be true if there is no spatial re- 
distribution of the energy released by Joule's effect inside the 
wire, but might be inaccurate if the wire height is not small 
compared to its width or if the wire thermal diffusion con- 
stant is very large compared to that of the substrate. In the 
extreme case where the wire temperature is uniform, simple 
calculations can be done, to provide an upper bound on the 
error. More precisely, for e ^ W and for times much larger 
than tw, we findQ j{x) 2^/{7rW^/l - {2x/W)^). Us- 
ing Eq. [3 we then find that the equilibrium temperature is 
lower than in the case of a uniform j{x) by an amount smaller 
than 5% for all the parameters explored in this article. Thus, 
we expect that the error in our predictions caused by our as- 
sumption on j{x) cannot be larger than 5%. Alternatively, 
the small wire temperature inhomogeneity predicted by our 
model shows that both boundary conditions give similar re- 
sults. 



III. COMPARISON WITH MEASUREMENTS 

We now come to the experimental measurements and their 
comparison with the model. The chip is based on an AIN sub- 
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strate of thickness 600(50) /im and of size 25 mmx35 mm. 
Gold wires of height 3 jj^m are deposited by evaporation on 
the substrate, on top of a 30 nm thick-titanium adhesion layer. 
We measure the wire heating after a constant current is turned 
on by monitoring the wire resistance and using Eq.|7l The cur- 
rent supplies we use allow a current rise time of about 10 fis. 
The wire resistance is deduced from the voltage drop across it. 
The contribution Rc of the connecting wires and the contact 
resistances has to be subtracted from the measured resistance. 
To compute we measure the circuit resistance at a cur- 
rent low enough to produce negligible heating and subtract 
the contribution of the microwire, computed using the nom- 
inal wire dimensions and the gold resistivity value at 300 K. 
We perform heating measurements on two different wires : 
a 200 /im wide, 20 mm long-wire and a 7 /im wide, 3 mm 
long-wire. 

To investigate the heating on long time scales we use the 
200 /im wide-wire, in which we run a current of about 5 A. 
Three main cases are considered. First, the chip is laid on 
paper to thermally insulate it from the copper block. Second, 
it is maintained on the copper block, either in air or in vacuum. 
Finally, the chip is glued on the copper block with Epotek 
H77. 

In the case where the substrate is laid on paper, we com- 
pare the data to the theoretical model, assuming that no heat 
escapes the substrate (Ic oo). As shown in Fig. [Till we find 
agreement within a few percent for times lower than 2 s, pro- 
vided the AIN thermal conductivity is set to A = 128 W/(Km). 
This value is close to the reported values for AIN II19|] al- 
though somewhat smaller. The discrepancy may compensate 
for imprecisions on the other parameters (AIN heat capacity, 
gold resistivity, wire size and substrate thickness) which are 
fixed to their nominal values. In all the following, we use this 
fitted value of A. With our substrate te — 7 ms, so that for the 
times considered in Fig. [Till one expects the heating process 
to be well explained by the IDx model of Eq. [3l Using this 
simple model and including the dependence of gold resistivity 
with temperature, we obtain the short-dashed curve shown in 
Fig-Hni It agrees within 5% with the more complete calcula- 
tion. The discrepancy, which is about a constant offset, is due 
to the early 2D regime. At times longer than 2.5 s, the data 
show an excess heating compared to the model. We attribute 
this effect to the finite size of the substrate : after t = 2 s 
the heat has diffused over a typical distance \^^Di =1.1 cm, 
larger than the distance from the wire to the substrate edge 
(1.0 cm). 

In the case where the substrate is maintained on copper, the 
relevant new parameter is the thermal contact resistance a be- 
tween the substrate and the copper block. We have measured 
a in the following way. First, we heat the substrate to a tem- 
perature Ti 500 K by running current in a wire and letting 
the substrate thermalize for a few seconds, while keeping it 
isolated from the copper with some paper. Then, we suddenly 
remove the paper so that the substrate is in thermal contact 
with the copper block, with an applied pressure of 2 x 10^ Pa. 
We monitor the thermal relaxation of the substrate by record- 
ing the resistance of one of the chip wires in which a small 
current of 10 mA is constantly flowing. A fit to the function 



T = To + (Ti - To)exp{-t/r) yields the relaxation time 
r = ace = 0.82 s, which gives the thermal contact resistance 
a = 5.77 X 10~^ m^KTW. To our knowledge no measurement 
of a was previoulsy reported for such a small contact pres- 
sure at a Cu-AIN interface. A much better coupling could be 
obtained for pressures higher than 10^ Pa ll24ll . but such high 
pressures do not seem realistic in the context of atom chips, 
because usually a large free surface on the chip is required. 
The value of a also depends on the surface quality, so we need 
to mention that the substrate we use has a backside roughness 
amplitude of about 100 nm. Importing in the model the mea- 
sured value of cr, we compute the heating of the wire after a 
5 A current- step, with no adjustable parameters. The result, 
shown as a dotted curve in Fig. Hill agrees with the measure- 
ment to better than 10%, which is a good validation for our 
model. 

In real experiments, unless the chip is itself one of the walls 
of the vacuum chamber 1I22I1 . the contact surface between the 
chip and the heat sink is in vacuum. We have thus also studied 
the case of a chip maintained on a copper block in vacuum, at 
a pressure lower than 10~^ mbar. We find that the thermal 
contact resistance avac between the substrate and the copper 
is much larger than in air, meaning that the thermal coupling 
in room conditions is actually caused by the air present in the 
voids between the substrate and the copper block. To measure 
c^vac, we monitor the cooling of the substrate with the same 
protocol as in air, with the difference that, since the cooling 
time is larger than the thermalization time of the substrate (a 
few seconds), thermal insulation of the substrate in the initial 
heating stage is not needed. The measured relaxation time is 
Tvac = 26 s, corresponding to a^^c = 1.8 x lO^^m^KAV. 
We also measure the relaxation time in vaccuum for a chip 
laid on paper and find 82 s, a value several times larger than 
Tvac- This confirms that conduction to the copper is the main 
cooling mechanism for a chip held on copper, even in vacuum 
where the thermal resistance dvac is high, and that black-body 
radiation is negligible. In practice, since the stationary time ts 
is identical to the relaxation time r, the coupling to the copper 
block has negligible effect on the wire heating for times much 
smaller than r. In particular, within the time of our measure- 
ment (4s), we observed a heating equal within 5% to the one 
observed when the chip is laid on paper. 

In order to mitigate the wire heating, especially if operating 
in vacuum, the previous studies show that it is highly desirable 
to fill the voids between the chip and its mount. That is why, in 
the last situation, we have glued the chip onto the copper block 
with a thin layer of Epotek H77, a glue of high heat conduc- 
tivity (according to the specifications, Xh77 = 0.66 W/(mK)) 
containing grains whose maximal radius is specified to 20 
fim. The experimental value of the steady state temperature 
(see Fig.HIB is obtained, within our model, for a thermal con- 
tact resistance a = 1.1 x 10~^ m^KTW, which corresponds 
to a glue layer of 73 /im, assuming a homogeneous glue layer 
and neglecting interface resistances. Importing this value into 
the model, we compute the expected wire heating, shown in 
Fig- nil as a dashed-dotted line. It matches the time evolution 
very well. 

We now turn to the measurements of the heating at early 
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FIG. 3: Wire heating at long times for a 200 /xm wide- wire running a 
current of 5 A. Solid lines give the experimental results for, from top 
to bottom, a chip laid on paper, a chip maintained on a copper block 
with a pressure of 2 x 10^ Pa and a chip glued on the copper block. 
The other lines are theoretical predictions (see text). 



times. With a 200 fim wide- wire, the early 2D heating regime 
is barely visible, because the condition <C e is not fulfilled. 
We thus use the 7 /im wide- wire. The measurements were 
taken for a substrate laid on paper and are shown on Fig. (H 
The prediction of the model with no free parameters, taking 
into account the 10 /iS current rising time, is shown as solid 
line. For times between 100 /iS and 5 ms, we clearly observe 
a logarithmic increase of the temperature, as expected from 
the 2D model, and the calculation agrees within 10% with 
the measurement. The dashed line is the result of the bare 
2D model of Eq. [2l including the temperature dependence of 
gold resistivity. Weff is adjusted to 0.15 W, so that the bare 
2D model and the more complete model agree within 1% at 
times smaller than 5 ms. At larger times, the model predicts 
a stronger heating due to the onset of the IDx regime, which 
is also visible in the experimental data. For times larger than 
50 ms, the model overestimates the temperature increase rate. 
This is very likely due to the finite length of the wire : indeed 
at t = 100 ms, the heat has spread in the substrate over a 
typical distance y/Dt = 2.3 mm that is no longer negligible 
compared to the wire length L = 3 mm. At longer times, the 
spread of the heat becomes two-dimensional in the plane of 
the substrate, which is more efficient for heat removal. 

Because it is experimentally relevant, it is worth investigat- 
ing briefly the effect of the finite wire length L on the sta- 
tionary temperature. As long as <C L, the equilibrium 
temperature is barely affected by the wire finite length. On 
the contrary, in the limit Is ^ L, the model of an infinite 
wire fails. In this case, for times t L'^/D, the wire tem- 
perature increases as AT2^^^(t) ^ {^L/47reX)ln{Dt/L'^) 
until the equilibrium temperature is reached. Equating the in- 
put power with the heat transferred to the heat sink, we find 
that the equilibrium temperature is about T^^^^ = T^^^ x 
(L/ttIs) ln{ls/L), a value much smaller than T^^^ (see sec- 
tion 2). Quantitative predictions are obtained using a two- 
dimensional version of Eq. [S] setting to and integrating 
over wave vectors in the x and y directions. In the experimen- 
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FIG. 4: Wire heating at long times for a 7 fim wide- wire running a 
current of 1.7 A. Lines are predictions of the full model (solid) and 
the bare 2D model of Eq. [2] (dashed). The crosses are experimental 
data. 



tal case considered in this paper, for a chip glued to the copper 
block and for our 3 mm long-wire (for which Ig = 0.96L), T^q 
is decreased by about 25% compared to its value for an infinite 
wire. 

The data from Fig. [4] can additionally be used to place an 
upper limit on the thermal contact resistance aw between the 
wire and the AIN substrate. In presence of a finite aw , 
the wire temperature increases by aw^/W on a time scale 
Tw = cTwCAuh, typically smaller than 1 fis, where Cau is 
the heat capacity of gold. This fast heating has been observed 
in Si-based atom chips ifTsll . Our data however are compat- 
ible with a vanishing aw, since they agree with our model 
which assumes aw = 0. More precisely, within the pre- 
cision of our model and measurements, we can confidently 
say that the excess heating due to a finite aw cannot be more 
than 10 K in the first 100 /iS, which yields the upper limit 
aw = 2.3 X 10~^Km^AV, corresponding to a conductance 
of 4.3 X 10^ W/(Km^). This contact resistance is a factor 6.6 
smaller than the values reported in Il8ll . 



IV. PRACTICAL CONSEQUENCES 

In order to minimize the wire temperature, it is of course de- 
sirable to minimize the thermal contact resistance a between 
the substrate and the reservoir. However, once a technology 
is chosen, so once a and Ic are fixed, it is possible to use our 
model to optimize the chip design and to compute the maxi- 
mal currents that can flow in each wire. 

As we have shown from our measurements, in good oper- 
ating conditions, a stationary state is expected to be reached 
within a few 100 ms. This means that typical experiments 
on atoms will fall into the stationary regime; for example, an 
evaporation to BEG typically takes a few seconds in our set- 
up. Therefore the stationary temperature ATgg is the relevant 
parameter to consider. In Fig. [51 a we plot ATgg for differ- 
ent wire widths and substrate thicknesses. The scalings used 
are justified by dimensional analysis : since heat conduction 
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inside the substrate is governed by linear equations, the equi- 
Hbrium temperature is given by 



a) 



ATeq = jf{W/lc,e/lc) 



(8) 



Note that the function / does not depend on c, which only en- 
ters into account for the time scales of the transient regimes. 
One sees in Fig.[5]a that ATeq is very well described by Eq.Hl 
in its domain of validity (1^ <C e <C Ic- Here Wq = 0.61iy 
has been obtained by fitting the calculation for = 10~^/c 
and 5 x 10"^ < e/lc < 0.3. The failure of Eq.|4]for the 
lowest curve of Fig. [5] a at small e/lc is due to the fact that for 
these parameters W is no longer small compare to Is. Then, 
the IDx regime barely exists and only the initial IDz regime 
is present. Roughly speaking, in the case e <C /c, the tem- 
perature distribution has a ID character, as can be seen in 
Fig. [5lb(B) which shows the temperature distribution inside 
the substrate for the parameters of Table 1 . In this regime, it 
is favorable to increase e, to "deconfine" the energy spread. 
On the other hand, for e ^ Ic, the IDx regime of diffusion 
no longer exists : the steady state is realized directly after the 
2D regime when the energy reaches the substrate's lower sur- 
face. ATeq is then of the order of \n{e/W)^/7TX, increasing 
with the substrate thickness (Fig.[5la). These arguments pre- 
dict an optimal substrate thickness of the order of Ic. This 
is confirmed by the calculations (Fig. [51 a), which show that 
ATeq reaches a minimum for the optimal substrate thickness 
e* :^ 0.31c, which is about independent of Vl^. In this opti- 
mal case, the three length scales Is, Ic and e are almost equal. 
In fact, e* is the smallest value of e that permits the suppres- 
sion of the IDx regime and therefore the temperature distri- 
bution inside the substrate has a two-dimensional character 
(Fig. [5lb(A)). Note that the wire heating minimum is rather 
broad (Fig. [51 a) : even with a substrate thickness as low as 
e*/10, the heating is increased by only about 50%. 

The second important quantity we can derive from our 
model is the maximal current that can flow in the wires. Using 
Eq.[71and <l> = pP/Wh, we find that Eq.[8lis a self-consistent 
equation, whose solution is ATeq = Po^ /{Wh)f /{I 
flcP/Wll) where Jq 
the critical current Ic - 



X^/ah/pQa. When / approaches 



Iq^W/ (Icf) a divergence occurs, 
signature of an instability. In Fig.[6l we show the critical cur- 
rent Ic, where the dependence on vW has been removed by 
dividing by ^JW/lc. Note that this maximal current corre- 
sponds to an established stationary state. For situations where 
currents are only required for a time shorter than ts, higher 
currents can be used. 



V. CONCLUSION 

To conclude, we have presented the first study on the ther- 
mal properties of AlN-based atom chips. The main result is 
that, as expected, the thermal behavior is more favorable than 
with Si chips. More precisely, our measurements are compat- 
ible with the absence of thermal contact resistance between 
the wire and the substrate. The heating of the wire is entirely 




b) 




FIG. 5: Optimum substrate thickness, a) Equilibrium wire tempera- 
ture versus substrate thickness for wire widths (from top to bottom) 
W/lc = 10"^, 10"^ 1.1 X 10-^ and 0.1. The dashed lines cor- 
respond to Eq. [31 with Weff = 0.61 VK. Heating is minimum for 
a substrate thickness of 0.3 Ic. b) Stationary temperature distribu- 
tion for the parameters of our chip glued on copper, for a substrate 
of thickness e=0.6 mm (B) and for a substrate of optimal thickness 
e* = 7.8 mm (A). The wire, 200 fim wide, runs a current of 5 A and 
its temperature increase is 24.3 K in case B and 15.8 K in case A. 




FIG. 6: Critical current divided by ^/W/lc and normalized to Iq = 

Xyah/poa, versus wire width for substrate thicknesses e/lc = 0.5 
(thick line, optimal case), 10 (dashed line) and 0.01 (thin line). The 
point represents the case studied in this paper : a 200 /xm wide, 3 //m 
high- wire on an AIN substrate of width e = 600 /xm, for the param- 
eters of Table 1. In this case, Ic = 16 A. 
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explained by heat diffusion inside the substrate and its absorp- 
tion by the heat reservoir holding the chip. We have devel- 
opped a model that accounts well for our experimental data 
for wires of different widths, for different couplings to the 
copper block and within the different heating regimes. The 
thermal coupling to the heat sink holding the chip is a crucial 
parameter : in particular we have shown that, when operat- 
ing in vacuum, it is recommended to glue the chip to the heat 
sink. Finally, the model is used to derive pratical learnings : 
the optimum substrate thickness is computed as well as the 
maximum current that can be run into the wires in the station- 
ary regime. 

We have treated a simplified case, but, in practical situa- 
tions, it is important to take other effects into account. First, 
as we already noticed, the finite length of the wire may reduce 
significantly the heating. The vicinity of wires may also af- 
fect the thermal behavior : several current-carrying wires sep- 
arated by distances on the order or smaller than the stationary 
length Is, will experience a stronger heating. In the temporal 
domain, the accumulation of heat over several experimental 



cycles will increase the heating if the dead time when currents 
are off is not long compared to ts . Finally, heat conduction 
through the copper block and its cooling may affect the heat- 
ing on long timescales. Active cooling of the copper block 
may be envisioned. 

Different technologies may improve the thermal properties 
of atom chips. An expensive solution would be to use dia- 
mond substrates because diamond, although electrically in- 
sulating, has a much higher thermal conductivity than AIN. 
Another key point is the realization of a better thermal con- 
tact between the substrate and the heat sink. If the backside 
of the chip is in air, using a thermal grease seems highly de- 
sirable. Working in vacuum, a process leading to a thinner 
glue layer may be developped and/or a different glue may be 
used. One could also consider soldering, the difficulty being 
to avoid damaging the chip. 
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